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This paper describes an algorithm for atmospheric state estimation based on a coupling between inertial navigation 
and flush air data-sensing pressure measurements. The navigation state is used in the atmospheric estimation algorithm 
along with the pressure measurements and a model of the surface pressure distribution to estimate the atmosphere using 
a nonlinear weighted least-squares algorithm. The approach uses a high-fidelity model of atmosphere stored in table- 
lookup form, along with simplified models propagated along the trajectory within the algorithm to aid the solution. 
Thus, the method is a reduced-order Kalman filter in which the inertial states are taken from the navigation solution and 
atmospheric states are estimated in the filter. The algorithm is applied to data from the Mars Science Laboratory entry, 
descent, and landing from August 2012. Reasonable estimates of the atmosphere are produced by the algorithm. The 
observability of winds along the trajectory are examined using an index based on the observability Gramian and the 
pressure measurement sensitivity matrix. The results indicate that bank reversals are responsible for adding 
information content. The algorithm is applied to the design of the pressure measurement system for the Mars 2020 
mission. A linear covariance analysis is performed to assess estimator performance. The results indicate that the new 
estimator produces more precise estimates of atmospheric states than existing algorithms. 


Nomenclature 


= backward smoothing gain 

= linearization of f with respect to x 

=  low-fidelity atmospheric model equations of 

motion 

= linearization of f with respect to u 

= gravitational acceleration, m/s? 

linearization of h with respect to x 

= pressure distribution model, Pa 

= identity matrix 

= linearization of h with respect to u 

= integer time index 

= integer time index of final pressure measure- 

ment 

= covariance of x after the measurement model 
update 

= pressure measurement vector, Pa 

= static pressure, Pa 

= process noise spectral density 

process noise covariance 

= pressure measurement covariance matrix 

= specific gas constant, J/(kg - K) 

= pressure measurement covariance matrix 

augmented with navigation uncertainty 
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S = prior covariance of x from low-fidelity model 

T = atmospheric temperature, K 

T = prior covariance of x from high-fidelity model 

t = time,s 

u = vehicle inertial state 

Uns Ves Ug = vehicle planet-relative north, east, and down 
velocity components, m/s 

W, = discrete-time observability Gramian 

Wnys Wes Wa = north, east, and down wind _ velocity 
components, m/s 

X11, X12, X29 = Van Loan integral submatrices 

x = atmospheric state vector 

At = time between pressure samples, s 

€ = pressure measurement error vector, Pa 

n = process noise input 

0, by = vehicle pitch, roll, and yaw attitude angles, rad 

P = density, kg/m? 

Omin = minimum singular value of the observability 
Gramian 

@ = state transition matrix 

Q = covariance of u 


I. Introduction 


ASA has developed an entry, descent, and landing (EDL) 

technology development roadmap [1] to guide investment 
strategies for increased EDL capabilities and robustness. One area of 
emphasis is on the development of precision landing capabilities 
achieved through improved environment/atmosphere characteriza- 
tion and EDL instrumentation for validation of engineering models 
and ground testing procedures. One approach that can be used to 
address these areas is the implementation of a flush air data-sensing 
(FADS) system, which uses an array of pressure ports installed in the 
vehicle forebody to measure the pressure distribution during entry. 
These pressure measurements can be processed to estimate the 
freestream aerodynamic state (such as flow angles, Mach number, 
and dynamic pressure), atmospheric conditions (density, pressure, 
and winds), and vehicle aerodynamics. These sensors can be used for 
postflight trajectory reconstruction and model validation, but they 
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also have the potential to be used to augment the onboard flight 
control system by providing estimates of density and winds if the data 
processing algorithms can be implemented in real-time. 

The traditional approach for estimating freestream atmospheric 
conditions from FADS systems is based on a multitiered approach in 
which first, the air data state (typically dynamic pressure, Mach 
number, angle of attack, and angle of sideslip) is estimated from the 
measured pressures, and then second, is combined with an inertial 
navigation system (INS) to estimate the freestream atmospheric 
conditions (density and winds). One example is from the space 
shuttle, in which air data from both the FADS system [2] and a 
deployable boom are combined with the INS to produce estimates of 
winds [3]. In this approach, the differences between the INS-derived 
flow angles (planet-relative) and the air data-derived flow angles 
(wind-relative) are attributed to winds, which results in a system of 
nonlinear algebraic equations. By neglecting the vertical wind 
component, the problem reduces to two equations (angle-of-attack 
and sideslip differences) and two unknowns (north and east wind 
components), which can be solved numerically using a Newton— 
Raphson method. The method has several advantages in that it is 
simple and yields a direct estimate of winds from in situ 
measurements based on the combination of inertial velocity, attitude, 
and wind-relative flow angles. Limitations of the approach are that 
the method is singular for flight-path angles of zero, and all 
instrumentation errors are mapped into the wind estimates rather 
being accounted for using an uncertainty model. 

Several approaches for tighter coupling between the FADS system 
and the INS have been proposed. One common technique is the use of 
complimentary filters [4,5], in which the attitude rates from the INS 
are used to smooth the FADS attitude estimate. The complementary 
filtering approach does not directly address estimates of the 
atmospheric conditions along the trajectory. Kalman filtering 
approaches have been proposed for wind and atmosphere estimation 
from the INS and air data sensors in [6-16]. Kasich and Cheng [6] 
developed an INS/FADS blended state estimator for computing 
freestream static pressure and aerodynamic flow angles, but the wind 
components were not estimated. A similar approach was proposed in 
[7] using a Kalman filter technique to estimate body-axis velocities 
and pressure altitude. Conversely, the studies of [8-12] developed 
Kalman filtering approaches suitable for real-time implementation, 
but they focused solely on estimating the atmospheric winds and did 
not address estimating the thermodynamic state of the atmosphere. 
Arain and Kendoul [13] developed an algorithm for estimating 
longitudinal winds for improved flight control. Karlgaard et al. [14— 
16] proposed full-state Kalman filters for estimating freestream 
winds, pressure, and density, but they were designed for postflight 
data processing rather than for real-time applications. 

In [16], a novel INS-aided FADS state estimation algorithm was 
developed in which the INS velocity magnitude was used in 
conjunction with estimates of static pressure and density to compute a 
Mach number pseudomeasurement. This approach provided 
enhanced FADS state estimates for flight conditions in the 
hypersonic regime, where estimating the Mach number from 
pressures alone was numerically problematic. The algorithm did not 
directly estimate freestream winds from the pressure measurement 
data; instead, wind components were estimated using the approach 
proposed in [3]. 

In this paper, the concept of the INS-aiding approach for the FADS 
state estimation algorithm is extended to include attitude information 
in addition to velocity. This approach enables the estimation of winds 
directly from the FADS pressure measurements by reformulating the 
state to include winds, density, and pressure. A reduced-order 
iterative Kalman filter method is implemented to blend information 
about the atmospheric state from the pressure measurements with 
prior knowledge of the atmosphere from models. These models can 
include tabulated high-fidelity models and simplified models for 
propagating the atmospheric state forward in time from one pressure 
measurement to the next. The algorithm is simple enough that real- 
time data processing is a possibility. A backward smoothing 
algorithm can also be used in postflight data-processing scenarios to 
provide improved estimates over the entire trajectory. 


The remainder of this paper is organized as follows. Section II 
describes the proposed atmospheric state estimator, including models 
and the filter/smoother algorithms. Section III applies this algorithm 
to the processing of the FADS data from the recent Mars Science 
Laboratory mission, and it compares results with previously 
published data from [16]. Section IV describes how this algorithm 
can be used to benefit the Mars 2020 mission, including optimization 
of the pressure port layout and pressure measurement accuracy 
requirements development. Concluding remarks are given in Sec. V. 
Note that the main interest in this paper is for entry vehicle air data 
processing; however, the algorithm is general enough that it could be 
used for wider applications such as high-speed aircraft. 


II. Atmospheric State Estimation 


The new approach for FADS state estimation proposed in this 
paper is to use the full state from the onboard navigation system 
(velocity and attitude) combined with the pressure data to directly 
estimate the atmospheric conditions. The estimation algorithm is 
aided by both high-fidelity atmospheric models that are tabulated vs 
altitude and simplified atmospheric models that are propagated along 
the trajectory within the algorithm. 


A. Pressure Modeling 


The processing of pressure measurement data to produce estimates 
of the freestream atmospheric conditions requires a mathematical 
model of the pressures at each port location as a function of the 
atmospheric conditions. This pressure model can be based on one of 
several different approaches, such as modified Newtonian flow [2], 
potential flow theory [17], or tabulated computational fluid dynamics 
(CFD) solutions [16]. In each case, these models can include 
correction factors that are used to match the pressure models to wind- 
tunnel test data. The pressure models are typically written as a 
function of angle of attack, angle of sideslip, Mach number, and 
dynamic pressure (or static pressure). In this paper, a new approach is 
taken and the form of the pressure model is recast to be a function of 
the inertial state of the vehicle (velocity and attitude) and the 
freestream atmospheric conditions (static pressure, density, and 
winds). Explicitly, the assumed model is of the form 


Px = A(Xp, Uy) + €y (1) 


where p, are the pressures at each FADS port at time k; x, = 
[P; Ps: Wn, We, Wal’ is the freestream atmospheric state; u,z = 
[U,, Ve. Ua, P, 9, w]' is the vehicle velocity and attitude state; and €, is 
the combined pressure transducer measurement error and pressure 
model error, which is assumed to be a zero mean with covariance R,. 
The rationale behind this assumed error model is that the pressure 
transducers are suitably calibrated from ground test data such that 
they are unbiased and the measurement covariance matrix is based on 
measured errors and quantifiable system uncertainties. 


B. Low-Fidelity Atmosphere Models 


A low-fidelity model of the change in atmospheric conditions 
along the trajectory can be derived from basic idealized relations such 
as the hydrostatic equation and the ideal gas law. Such simplified 
relationships are suitable for implementation in the algorithm for 
propagating the atmospheric state estimate forward between pressure 
measurements, which are assumed to occur at a reasonably high rate 
(several samples per second) along the trajectory. Since the simplified 
model involves idealized approximations, uncertainties in the model 
can be accounted for with process noise. 

A model for the rate of change in static pressure can be found by 
rewriting the hydrostatic equation as the time derivative of pressure 
along a given trajectory, namely, 


Ps = pP8vq (2) 


Similarly, a model for the rate of change in density along the 
trajectory can be derived from the ideal gas law, with the assumption 
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that the atmosphere is locally isothermal (T = 0) between FADS 
pressure samples. The equation is of the form 


Ps = PsP = gvup” 
RT Ds Ps 


p= (3) 


A commonly used model for the rate of change in atmospheric 
wind is a random walk model [12,14-16,18,19] where the 


deterministic portion of the model is simply w, = w. = Wy = 0. 
Thus, the low-fidelity model can be written in the form 


x= f(x.u)+n (4) 


where 


gvap’/Ps 
P8vVa 
f(x, u) = 0 (5) 
0 
0 


The quantity 7 in Eq. (4) is a process uncertainty term that is 
assumed to be zero mean with spectral density Q. This term accounts 
for uncertainties in the simplified model, and it is a trajectory- 
dependent tuning parameter. 

The continuous model in Eq. (4) can be transformed to a discrete 
model of the form 


Xeq1 = Xe + L(y, Uj )At (6) 


which is suitable for propagation between pressure measurements. 
The uncertainties in the model can be propagated between pressure 
measurements using the relation 


Siat = ®, 5,0) + GQ,G] + 0; (7) 


where ®, is the state transition matrix, and 0; is the discrete-time 
process noise covariance. These quantities can be jointly calculated 
from the van Loan matrix integral [20], given by 


—F, OQ —[Xn Xo] _[X, &;'@ 
ero(| 0 m |e) =| 0 2] =| 0 | » 


which leads to the results ®, = X}, and 0, = ®,X,,. Assuming a 
reasonably small integration time step, these quantities can be 
approximated by ®, ~ I + F,,At and Q; = Q,At. 


C. High-Fidelity Atmosphere Models 


High-fidelity atmosphere models can also be incorporated into the 
state estimate as prior information. These atmosphere models are 
generally computationally intensive enough that implementation 
within the algorithm is infeasible. Instead, the high-fidelity 
atmosphere model data can be incorporated using table lookups 
where the atmospheric conditions and uncertainties are tabulated as a 
function of altitude along some nominal trajectory. The model of this 
form produces an estimate of the atmospheric conditions x,, along 
with an associated error covariance matrix T;. 


D. Data Fusion Algorithm 


The atmospheric state estimate can be determined from a fusion of 
the available data sources, including the FADS pressure 
measurements, as well as information from the low- and high- 
fidelity models. The algorithm is in the form of a nonlinear weighted 
least-squares solution with the model data incorporated as prior 
estimates. 

The pressure measurement model can be approximated by means 
of the truncated series expansion 


Dy & W(X_, Up) + Ay (xy — Xp) + Ex (9) 


where x; is some reference state, H;, is the measurement sensitivity 
matrix, and €, is an error state that includes pressure measurement 
and model uncertainties in addition to uncertainties in the navigation 
state u,. The covariance of €; is given by R, =R,4+ JQ J]. 

The atmospheric state estimation problem can be reduced to a 
linear regression problem of the form y, = Hyx, + €,, where 
VY. = Py —h(xy, u,) + A,x,. By virtue of the Gauss—Markov 
theorem, the best linear unbiased estimate of x,, denoted by xX, is the 
weighted least-squares solution [21], where the atmospheric model 
data are treated as prior observations 


&, = (ATR! Ag + Sp! + Tz!) (AT Re ye + Se + Tz! xn) 
(10) 


Since the relationship between the states and the measurement is 
nonlinear, the estimation scheme can be iterated until convergence by 
successively replacing x, by X,. The state estimate error covariance 
matrix P; of the converged solution can then be computed from 


P, = (ATR, A, + So +7;')"! (11) 


After incorporating the pressure measurement data and prior 
information from the atmospheric models, the best estimate of the 
atmospheric state and its covariance can be propagated forward in 
time to the next measurement sample using the low-fidelity model 
equations of motions given in Eqs. (6) and (7). Note that this 
framework is essentially an iterative extended Kalman filter, with the 
additional incorporation of prior data from the high-fidelity model 
tables. In a sense, the proposed technique is an optimal blending of 
FADS pressure measurements and a virtual air data system [22] that 
uses only forecast atmosphere data. 


E. Backward Smoothing 


The estimate of the atmospheric state at the end of the process has 
incorporated all available data, and thus represents the best possible 
estimate of the atmosphere. In a postflight processing scenario, this 
best estimate can be smoothed backward to the initial time in order to 
map this information over the entire trajectory. This process is known 
as fixed-interval smoothing. One well-known solution of the fixed 
interval smoothing problem is the Rauch—Tung-—Striebel smoother 
[23]. The backward recursion formulas are given by 


Cy = Py P] Pepi (12) 
Ke = Xen + Ceri — Xie — LEies MAM (13) 
Puy = Pye + CP esi — Presid] (14) 


The backward smoothing procedure begins with the final state 
estimate given all available measurements Ey n and its covariance 
Py, and it propagates from k = N backward to k = | using the 
relations given previously. 


F. Aerodynamic State Transformations 


The atmospheric state (winds, pressure, and density) are outputs of 
the proposed FADS data processing algorithm (from either the 
forward filter for possible real-time applications or the backward 
smoother for postflight processing). The atmospheric state can 
readily be combined with the INS state solution (velocity and 
attitude) to produce estimates of aerodynamic states, including angle 
of attack, sideslip, Mach number, dynamic pressure, and (when 
combined with mass properties) vehicle aerodynamic coefficients. 
Uncertainties can be mapped from the atmospheric and INS states 
into the aerodynamic states through a linear covariance analysis. The 
equations of the transformation from atmospheric and INS states to 
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aerodynamic states are readily available in various sources such as 
[24] and are not repeated here. 


II. Application to Mars Science Laboratory 


On 5 August 2012, the Mars Science Laboratory (MSL) entry 
vehicle successfully entered the atmosphere of Mars and landed the 
Curiosity rover safely on the surface of the planet in the Gale crater. 
The MSL entry vehicle was comprised of a 70 deg sphere-cone 
heatshield and backshell consisting of a stack of three truncated 
cones. The forebody was similar to the heatshield geometry 
developed for the Viking Mars landers. A phenolic-impregnated 
carbon ablator (PICA) was used for the thermal protection system 
material. The backshell configuration was also similar to Viking, 
with a third cone section added to accommodate the parachute 
volume. The MSL vehicle as-built outer mold line is shown in Fig. la. 
During most of the entry, the capsule used a radial center of mass 
offset to fly at an angle of attack (approximately 16 deg at hypersonic 
conditions). This attitude produced lift to fly a guided entry profile, 
reducing the landing footprint to a much smaller size than any 
previous Mars mission. To fly the guided entry, the vehicle carried 
four pairs of reaction control system jets to perform maneuvers and 
damp rates. The four pairs of jets could be fired rapidly in different 
combinations to provide control torques about roll, pitch, yaw, or any 
other axis by modulating the pulses of the jet. 

The MSL carried with it an instrumentation package designed to 
measure the aerodynamic and aerothermal environments during 
atmospheric entry. This instrumentation package was known as the 
MSL entry, descent, and landing instrumentation (MEDLI) [25], 
which consisted of three major subsystems: the Mars Entry 
Atmospheric Data System (MEADS), the MEDLI integrated sensor 
plugs (MISPs), and the sensor support electronics (SSEs). The 
MEADS consisted of seven pressure transducers connected to flush 
orifices in the heat shield to measure pressures across the vehicle 
forebody. The MISP devices were a system of seven thermocouple 
and recession sensors that provided aerothermal measurements of the 
heat shield performance. The SSE provided power to the sensors, 
conditioned their signals, and transmitted the data to storage on the 
Curiosity rover. The MEDLI sensors provided measurements that 
were used for trajectory reconstruction and engineering validation of 
aerodynamic, atmospheric, and thermal protection system models in 
addition to Earth-based systems testing procedures. The MEDLI data 
and their usage for reconstructing the aerodynamic and aerothermal 
performance of the MSL entry vehicle were described in [16,26-30]. 


A. Flight Reconstruction 


This section describes the results of applying the new full state 
INS-aiding approach to the MEADS atmospheric reconstruction. 


The high-fidelity atmosphere models are based on combined 
mesoscale models, which include the Mars Mesoscale Model, 
Version 5 [31]; and the Mars Regional Atmosphere Model System 
[32]. The mesoscale atmosphere models provided estimates of the 
mean and variability of atmospheric values from the surface to an 
approximately 50 km altitude. Further details of how the models were 
used for MSL can be found in [33]. Uncertainties in the INS solution 
were based on linear covariance propagation of the MSL navigation 
filter algorithm [34]. 

Figures 2-4 show a comparison of the basic (unaided) FADS 
algorithm and the velocity-based INS-aided FADS algorithm from 
[16]. Both algorithms produce good estimates of dynamic pressure 
but, due to the inherent weak observability of static pressure from 
surface pressure measurements alone [16], the unaided FADS 
algorithm produces wildly varying estimates of the Mach number. As 
a consequence, there is some angle-of-attack error buildup in the 
unaided FADS algorithm because the algorithm assumes an 
erroneous Mach number, and thus looks up pressures in the CFD 
database at an incorrect flight condition. The velocity-based INS- 
aiding algorithm stabilizes the estimates of Mach number and static 
pressure. The new proposed tightly coupled INS/FADS algorithm 
produces estimates that are similar to the velocity-based INS-aiding 
algorithm. 

Figure 4 shows a comparison of the density and static pressure 
estimates from the three algorithms. Although the density estimates 
are comparable between the three algorithms, the static pressure 
estimate from the unaided FADS algorithm is highly erratic and 
completely unreasonable due to the inherent weak observability of 
this method. Both the aiding algorithms produce reasonable 
estimates of static pressure. Some slight differences are apparent 
between the two aiding algorithms in the estimate of atmospheric 
density. Here, differences arise due to the treatment of winds in the 
two algorithms. The velocity-aiding algorithm assumes zero winds, 
and it computes density directly from dynamic pressure and INS 
relative velocity. Thus, nonzero winds are a possible error source in 
the algorithm that would result in a biased estimate of density. The 
tightly coupled algorithm produces a slightly different estimate due to 
the fact that winds are estimated from the pressure measurements, and 
thus the estimated density is compensated for winds. The difference 
between these density profiles is consistent with a tailwind, which 
would have the effect of reducing the apparent density if winds are 
otherwise assumed to be zero. 

The reconstructed atmospheric winds are shown in Fig. 5. These 
figures show comparisons between the winds derived from the 
velocity aiding solution and the new tightly coupled INS/FADS 
algorithm. Recall that the winds computed from the velocity-aided 
solution were obtained by the method proposed in [3], wherein the 
wind estimates were produced by reconciling differences between 
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Fig. 1 Vehicle geometry (BIP, backshell inferface plate; SC, spacecraft coordinates; BS, backshell; HS, heatshield; Sep., separation; IMU, inertial 


measurement unit; SR, spherical radius). 
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Fig. 4 Atmosphere. 


a) Density 


respectively. These estimates are reasonable given the atmospheric 
uncertainties [28], and moreover support circumstantial evidence 


FADS (wind-relative) and INS (planet-relative) angles of attack and 


for winds based on the vehicle entry guidance response [35] and 


sideslip. This algorithm produces wind estimate results that vary 


wildly over the trajectory and, for the most part 


trajectory [36]. The wind estimates are also consistent with 


, produces 


unreasonable results. The algorithm is better behaved at lower Mach 
numbers where winds are more observable because they are a larger 


portion of the total air-relative velocity. 


postflight full-state Kalman filter/smoother results from [16]. Note 


that, following the entry balance mass (EBM) jettison event, the 


pressure data become corrupted by structural vibrations due to 


The new tightly coupled INS/FADS algorithm produces 


reasonable wind estimates over the entire entry trajectory. These 
wind estimates indicate a north wind component of approximately 
—10 m/s and an east wind component of approximately 20 m/s 
over much of the entry trajectory. For the MSL entry trajectory, the 
north and east winds are essentially pure cross and tail winds, 


pyrotechnic shocks. These shocks cause the transducer diaphragms 


to vibrate, which increases the noise in the pressure data. This 
increased noise level shifts the filter weighting such that the wind 
estimates revert to the tabulated atmosphere model following the 


EBM jettison event. 
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Fig.5 Winds. 


B. Observability Analysis 

The MSL entry vehicle flew a guided entry, using bank reversals to 
control downrange position during hypersonic flight [35]. A time 
history of the bank profile and resulting vehicle yaw attitude is shown 
in Fig. 6a. These bank reversals are designed for the purpose of EDL 
guidance, but they have the secondary benefit of increasing the 
observability of winds along the trajectory. The observability of the 
wind states along the trajectory can be assessed via the discrete-time 
observability Gramian, defined as [21] 
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One scalar measure of the observability of the system is the 
minimum singular value of the observability Gramian [37]. If the 
minimum singular value is small, the system is difficult to observe. 
Likewise, large values of the minimum singular value indicate the 
opposite. Other observability metrics include the rank, determinant, 
eigenvalues, maximum singular value, trace, and condition number 
of the observability Gramian [37]. The minimum singular value of the 
observability Gramian o,,;, and its time rate of change along the MSL 
entry trajectory are shown in Fig. 6b. By examining the rate of change 
of Opin, It is apparent that the minimum singular value starts to 
increase just after the time of the bank reversals. The implication is 
that the bank reversals have the effect of adding information content 
to the system. 
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Fig. 6 Bank profile and observability. 
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b) Sensitivity to pressure 


Fig. 7 Pressure measurement sensitivities to atmosphere states. 
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angle such that winds can be determined easily from the sideslip 
pressure port measurements. Head/tail winds are more difficult for 
the filter to resolve because the centerline ports are strong functions of 
both the angle of attack and dynamic pressure, which serves to couple 
the head/tail wind with density as well as the downward wind 

The observability of winds due to the bank reversals has several 
consequences. First, a purely ballistic entry with no guidance may not 
be able to resolve winds along the trajectory. Second, the sensitivity 
of winds to bank reversals opens up the possibility of designing 
maneuvers to maximize the observability of winds along the 


trajectory. A crosswind corresponds to a nonzero inertial sideslip 
component. 


northerly wind corresponds 
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correlation to the bank profile. The implication is that the bank 


maneuvers add information content about the atmospheric winds. It 
can also be seen that the pressures are more sensitive to the northerly 
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The source of the information content can be gleaned by 
examining the pressure measurement sensitivities with respect to the 
states along the trajectory. Figure 7 shows the pressure measurement 
sensitivity with respect to density and static pressure. These plots 
show that the pressure measurements are highly sensitive to density, 
which is to be expected. The sensitivity with respect to density does 
not show any correlation to bank maneuvers. The static pressure 
sensitivity shows some correlation to the bank profile. 

The pressure measurement sensitivities with respect to the wind 
components are shown in Fig. 8. These sensitivities show a strong 
component of the wind. As stated in the previous section, the 
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Fig. 8 Pressure measurement sensitivities to wind states. 
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Fig. 10 Port arrangement. 


estimation has been proposed in [38], for example. Such a process 
would be akin to designing system identification maneuvers for 
estimating aircraft aerodynamic parameters. 


IV. Implications for the Mars 2020 Mission 


Another MEDLI-like system of instruments is planned to be flown 
on the Mars 2020 mission. This instrumentation system, known as 
MEDLI2 [39], will acquire FADS pressure data to be used for the 
reconstruction of atmospheric states and vehicle aerodynamics 
during entry. The focus of the pressure system on MEDLI23 is geared 
toward estimating aerodynamics in the supersonic regime of 
flight, where some questions remain regarding the aerodynamic 
reconstruction of MSL [29]. To this end, the forebody pressure 
system will carry one transducer with a full-scale range of 35 kPa (the 
same as MSL transducers) to measure stagnation pressure over the 
entire entry trajectory (which in tum yields estimates of dynamic 
pressure and density) and six transducers with a full-scale range of 
7 kPa to more accurately measure the atmosphere and aerodynamics 
in the supersonic regime of flight (roughly Mach 6 and below). In 
addition, one transducer will be installed on the backshell to measure 
the base pressure and its contribution to drag. Since the focus of this 
instrumentation is on supersonic measurements, estimates of winds 
are far more critical to interpreting the reconstructed aerodynamics. 


The diagram showing the MEDLI2 pressure port locations on the vehicle 
forebody was created by John Van Norman. 


Thus, the algorithm described in this paper is expected to prove useful 
in the postflight data reconstruction effort. 

The following two sections describe the design of the pressure port 
layout to optimize atmospheric state observability and the linear 
covariance analysis of the estimator performance for a given pressure 
measurement uncertainty. The Mars 2020 reference trajectory on 
which these analyses are based in shown in Fig. 9. 


A. Port Placement Optimization 


The reformulated wind estimation algorithm described in this paper 
enables a straightforward method to optimize the FADS port layout in 
order to provide the best state estimates for a given trajectory. The 
optimization can be achieved by solving for pressure port locations 
that maximize the observability (or, equivalently, minimize the 
unobservability) of the atmospheric state based on the observability 
Gramian. Sensor placement optimization for observability was first 
suggested in [40], and it has subsequently been studied for a wide range 
of applications such as structural vibration [41], chemical reactors [42], 
airfoil wake estimation [43], and Mars entry navigation [44]. 

Another approach, suggested in [45], involves sensor placement 
optimization that minimizes the trace of a weighted state covariance 
matrix. This approach has been used in the past for FADS pressure 
port location optimization in [46-48] (minimum root-mean- 
square errors were used in [46)). This approach has the drawback 
that the optimal sensor locations become dependent on modeling 
assumptions such as the process and measurement noise covariance 
and the choice of the weighting matrix. 
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Fig. 11 Pressure measurement sensitivities to dynamic pressure, density, and flow angles. 
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Fig. 12 Pressure measurement sensitivities to wind states. 
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In this paper, the pressure port locations are chosen to directly 
maximize the observability of the system by minimizing the 
reciprocal of the minimum singular value of the observability 
Gramian. In this method, the optimal port locations are functions only 
of the reference trajectory and the CFD pressure distribution model. 
A sequential quadratic programming method is used to solve the 
optimization problem. Several constraints are imposed on the 
optimization, namely, that the ports are constrained to lie within a 
boundary in order to not be placed near the shoulder, away from 
leeside turbulent conditions, and at least 76 mm (3 in.) away from 
PICA tile seams. Ports 1, 2, and 5 are constrained to lie at the 
stagnation regions before and after the entry ballast mass ejections, 
which occur shortly before parachute deployment and involve a 
change to the trim angle of attack. The remaining ports (3/4 and 6/7) 
are constrained to be in symmetric pairs across the vehicle pitch plane 
of symmetry. 

A candidate port arrangement for the Mars 2020 mission is shown 
in Fig. 10a. The port pressures along the trajectory in the supersonic 
regime are shown in Fig. 10b. Note that the supersonic pressure 
measurements are saturated above 7 kPa. 

Pressure measurement sensitivities with respect to the dynamic 
pressure, density, and aerodynamic flow angles are shown in Fig. 11. 
These results indicate that ports 1 and 2 provide the most information 
about the freestream dynamic pressure, followed closely by ports 3 
and 4. This result is expected due to the proximity of the ports to the 
forebody stagnation region. The same trend is evident in the pressure 
measurement sensitivity to freestream density. Ports 5, 6, and 7 
provide the most information about the angle of attack. The angle of 
sideslip is primarily sensitive to ports 6/7 and 3/4, which is also 
expected because these ports lie in symmetric pairs off the vehicle 
pitch plane of symmetry. 

Pressure measurement sensitivities to the wind states are shown in 
Fig. 12. Changes in the measurement sensitivity along the trajectory 
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show some relationship to the bank angle, as expected based on the 
discussion of wind observability described in Sec. III.B. Note that, for 
this particular reference trajectory, north and east winds do not 
correspond to cross and head winds, respectively, as was the case for 
the MSL. 


B. Linear Covariance Analysis 


This section describes a linear covariance analysis of the estimator 
performance. The uncertainties are propagated along a nominal 
trajectory, assuming pressure measurement accuracies of 1% of 
reading, and pressure port locations shown in Fig. 10a. CFD 
uncertainties and INS uncertainties are ignored for this analysis, 
simply to show a comparison between the velocity-aided algorithm 
developed in [16] and the new tightly coupled algorithm developed in 
this paper. 

Figure 13 shows the linear covariance analysis results for estimates 
of dynamic pressure and Mach number. The results indicate that the 
dynamic pressure estimates corresponding to the new tightly coupled 
algorithm are roughly the same precision as the velocity-aided 
algorithm at high Mach numbers, and the new algorithm becomes 
more precise than the velocity-aided algorithm as the Mach number 
decreases. The Mach number reconstruction corresponding to the 
tightly coupled algorithm is superior to the velocity-aided algorithm 
across the range of Mach numbers. This result is due to both an 
enhanced estimate of the speed of sound as well as an improved 
estimate of winds (which in turn produces an improved estimate of 
wind-relative velocity), which is shown in the following figures. 

A comparison of results for the aerodynamic flow angles is shown 
in Fig. 14. The results are fairly comparable, although the tightly 
coupled algorithm produces a slightly improved result over that of 
the velocity-aided algorithm. This improvement is due to the 
incorporation of additional information in the form of the INS attitude 
estimates. 
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Fig. 14 Linear covariance analysis results: aerodynamic flow angles. 
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Fig. 16 Linear covariance analysis results: winds. 


Uncertainty analysis results for the atmosphere states are shown in 
Fig. 15. These results indicate that the tightly coupled algorithm is 
able to produce improved estimates of the atmospheric states. The 
density estimate is improved because of improved estimates of the 
wind state. The velocity-aided pressure reconstruction uncertainty is 
a highly conservative estimate of the true uncertainty because this 
pressure estimate arises by integrating the hydrostatic equation, 
treating each density sample as uncorrelated with previous estimates. 
In reality, these estimates are expected to have some level of 
correlation, which would have the effect of reducing the pressure 
estimate uncertainty. 

Finally, wind estimate uncertainties are shown in Fig. 16. The 
tightly coupled algorithm produces improved estimates of winds by 
incorporating INS attitude measurements directly into the estimation 
process. 


V. Conclusions 


A new inertial navigation system aiding approach is developed for 
flush air data-sensing systems. The approach uses a full-state aiding 
(velocity and attitude) opposed to past approaches that used only 
velocity magnitude. The advantage of the new aiding approach is that 
winds can directly be estimated from the measured pressures, rather 
than being inferred from discrepancies between inertial angle of 
attack and wind-relative angle of attack that can sometimes lead to 
singularities or other issues with observability of winds. The 
algorithm produces a reasonable estimate of the atmospheric 
conditions based on flight data from the Mars Science Laboratory 
entry vehicle that was acquired in August 2012. In particular, the 
estimates of winds are consistent with circumstantial evidence based 
on vehicle dynamics and with full-order Kalman filter/smoother 
results. Linear covariance analysis results indicate that the new 


algorithm produces enhanced atmospheric state estimates when 
compared to the existing state-of-the-art air data processing 
algorithm. The algorithm is reasonably straightforward such that 
real-time implementation is possible, although computational 
complexity/feasibility would need to be assessed for each particular 
application. 
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